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ABSTRACT 



For studying the thermodynamic properties of systems using statistical mechanics we 
propose an ensemble that lies in between the familiar canonical and microcanonical 
ensembles. From a comparative study of these ensembles we conclude that all these 
ensembles may not yield the same results even in the thermodynamic limit except at 
high temperatures. An investigation of the coupling between systems suggests that the 
state of thermodynamic equilibrium is a special case of statistical equilibrium. As a 
byproduct of this analysis we have obtained a general form for probability density 
function in an interval. 



In studying the thermodynamic properties of systems using statistical mechanics, we 
restrict ourselves to a constant energy surface since we know that the energy of the 
system under consideration is constant. The basic assumption of ergodicity which 
helps us replace the time average of observables by the corresponding phase averages 
leads us naturally to a constant density on the energy surface given by 



The ensemble of systems represented by this equation is called the microcanonical 
ensemble in which all the systems have the same energy E. For most physical 
systems the use of microcanonical ensemble turns out to be formidable. Besides this 
Gibbs 1 proposed the canonical ensemble. Starting from the postulates that 1) for 
statistical ensemble density p has to be a stationary state of Liouville's equation 
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I. INTRODUCTION 



p = a constant if (E-AE)<E<(E + AE) and 



otherwise 



(i) 




(2) 



and that 2) 



(3) 



where H is the Hamiltonian and rfFis the volume element in phase space he 
proposed 



p = cxp\-a- PH\ 



(4) 



as "the most simple form conceivable". Exp {-a} is the normalisation constant in the 
above equation. He went ahead to identify /3 as l/(kT), where k is the Boltzman 
constant and T is the temperature, from the structure of a Pfaff differential which is 
analogous to that of dF in thermodynamics. The ensemble of systems represented by 
Eq. (4) is called the canonical ensemble of Gibbs. It was only many years later that 
Jaynes 2 showed that these ensembles can be derived from information theory (3 ' 4) . 
Starting from this theory we embark upon the intermediate ensemble given by, 

p = exp[-a -f3H- (y/N)H 2 ] (5) 

where jfN is the LaGrange multiplier corresponding to an additional constraint that 
has been introduced as a better approximation to the microcanonical ensemble, the 
subject and purport of this exercise. In section II we outline the method of 
information theory and derive all possible forms of density thereof. We also obtain 
general forms of probability density function 5 and compare them with those obtained 
from information theory. Since our aim is to study the thermodynamic properties of 
systems we identify, in section III the coefficients /3 and yfN of Eq. (5) in terms of 
thermodynamic quantities and at once realise the physical reasoning for introducing 
higher order terms of H in p. In section IV we derive the analogous thermodynamic 
quantities. From a comparative study of the mean energy and fluctuations in all 
these ensembles we conclude in section V that canonical and microcanonical 
ensembles may not yield the same results even in the thermodynamic limit except at 
high temperatures. The coupling between systems described by various ensembles is 
investigated in section VI and the conditions for statistical and thermodynamic 
equilibria discussed. We summarise the results in section VII. 

II. INFORMATION THEORY AND DENSITY FUNCTIONS 

The method of information theory is but the method of finding the most probable or 
least biased density function involving LaGrange' s multipliers. To find 5 the most 
probable distribution p subject to a given set of linearly independent constraints that 
fixes all the moments 

< X i >= C. (6) 

where kr s are constants and i takes all nonnegative integer values, assign a statistical 
entropy to p defined by 

S P = -(k In p) (7) 

In Eq. (6) the brackets denote the averaging done with respect to the density function 
and k in Eq.(7) is introduced for thermodynamic reasons. The least biased form of p 
is obtained by minimising information, which is nothing but negentropy, subject to 
the complete set of linearly independent constraints given by Eq. (6). For any small 
variation in p we have 



5\dT p [-k In p -kz™ a x ] = (8) 

which yields 

p (x) dx=exp [— iji Cfcx' ] dx (9) 

where C^-'s are LaGrange multipliers to be determined from Eq. (6). If the constraints 
are only n in number then this method gives the most probable least biassed form of 
density whose first n moments are fixed. For instance if we do not know anything 
about the system except that the particle number and the energy are fixed we get on 
the constant energy surface 

p =exp[-a] (io) 

which is the microcanonical ensemble. Exp( oi) gives the possible number of states 
with a given energy. If on the other hand we want to define the ensemble in such a 
way that the average value of the energy of the ensemble alone is fixed we obtain 



p = cxp\-a- /3H\ (id 

which is the canonical ensemble and 

p = exp[-a-/3H-(y/N)H 2 ] (12) 

would correspond to that ensemble obtained by conserving (H) and (H ). This 
ensemble is a subset of Gibb's ensemble and it encompasses the microcanonical 
ensemble. By conserving higher and higher powers of H we would approach 
microcanonical ensemble. Hence 

p = exp[-E- ja#'] (is) 

represents the microcanonical ensemble density and representing as it does the entire 
class of physical systems it must be nothing but a general form of probability density 
function in an appropriate interval and by truncating it any intermediate ensemble 
could be realised. 

For a probability density function of one variable Eq. (9) can be derived in yet 
another way 5 . Let p(x)dx be a continuous probability density function in [a,b] that 
has no real or imaginary zeros or singularities. From Taylor expansion of lnp(x) 
about Xo£[a,b] we get 

In p(x) = In p(xo) + £~ (x-x ) 1 1 / i ![-^ln p(x)] ( * = «> (14) 

dx' 

and from the above 

p(x)dx = exp[-Y J °° =Q adc']dx (15) 

which is of the same form as Eq.(9) If p(x)dx has real or imaginary zeros or 
singularities at a and/or b then the right hand side of Eq.(14) is not convergent and if 



we represent the zeros and singularities by Z(x) and P(x) respectively, Taylor 
expansion of ln{p/(Z(x)P(x))} still converges yielding 



p(x)dx = Z(x) P(x) exp [-Y,Zo aoc (16) 

Comparing Eq.(15) and Eq.(16) we realise that if information associated with the 
probability is redefined 5 as 

/ = (ln{ PW }) (i7) 
l Z(x)P(xy 

then the results obtained from information theory and Taylor expansion tally. If 
p(x)dx has a finite number of real or imaginary zeros or singularities or a finite 
number of finite discontinuities in the interior of the interval we only have to divide 
the given interval into subintervals to follow the above procedure. 

This can also be extended 5 to the case of many variables. For example a 
probability density function of two variables arising from 

(x l y jX j = dj (18) 

yields 

p(x,y)dxdy = exp[-£~ J)=(00) a^y ]dxdy (19) 

where (x,y)eDaR 2 and p does not vanish anywhere in D. Coefficients ay are derived 
from 

-d 

— {InpO)} ) = Cij (20) 
octa 

Eq. (20) also can be derived 5 from Taylor expansion. Vide appendix. 

III.IDENTIFICATION OF LAGRANGE MULTIPLIERS 

Having derived the formal expression for ensemble densities we now proceed to 
identify the LaGrange multipliers in terms of thermodynamic quantities. The 
transition from statistical mechanics to thermodynamics is effected in one step where 
the statistical entropy of the microcanonical ensemble is identified as the 
thermodynamic entropy. From Eq. (7) we get 

S P = k In g(E) (2i) 

where is g(E) is the possible number of states with a given energy E. From Eq.(21) 
the LaGrange multipliers occurring in Eq.(12) and Eq.(13) can be easily evaluated. It 



ensures, as will be shown below, that all the terms in the exponent in Eq. (13) are 
linear in N, the particle number. With this objective let us apply the microcanonical 
ensemble to the system and the rest of the system. Let Es and Er and Ns and Nr 
and Vs and Vr denote the energy, particle number and volume of the two respectively 
such that Es+E R +Ei nt = E where Ei„ t is the interaction energy much smaller than E& 
The entropy of the rest of the system is 

Sr(Er) = S R (E) + ^J-E s y li\ [d* (22) 

by Taylor expansion and 



ds* 


1 


d Eu 


~ kT 


d 2 S* 


-1 


dE R 2 


kT 2 Cv R 



where the temperature T and the specific heat Cy refer to the rest of the system. Since 
the system and the rest of the system are decoupled- Ei„ t «Es , the ensemble density 
for the system obeys the following equation. 



P oc g R (E R ) (25) 

If gR(ER) is large enough so that Er can be treated as a continuous variable we obtain 
from Eq. (21), Eq. (22), Eq. (23) and Eq. (24) 

p =exp[-Es/kT -Es 2 /(2kT 2 Cv) -....] (26) 

In the above equation whereas the first exponent varies as Ns the second varies as 

Ns /Nr. In the limit when (Ns/Nr)—X) the second and higher order terms in Es vanish 
thus yielding Gibbs ensemble. This limit defines the concept of infinite reservoir. An 
infinite reservoir is the requisite of an experimentalist. When he measures the 
temperature of a system he should measure it in such a way that the system is itself is 
not in the least disturbed. Needless to say that it need not be and should not be 
introduced in the theoretical definition of equilibrium. For measuring temperature 
experimentally we do need it. It must be added that Gibbs 1 brings in the concept of 
heat bath only while interpreting exp(-/3H) and not for defining the equilibrium state. 
Thus when we relax the condition (Ns/Nr)—>0 the higher order terms in Es would no 
longer be negligible. 

How do we decide what (Ns/Nr) should be? When we say that the system is in 
equilibrium at a temperature T, then every macroscopic subsystem is in equilibrium 
with every other equally large macroscopic subsystem at the same temperature T 
which would also ensure equilibrium with infinite heat bath. But the latter does not 
imply the former. Hence we define the equilibrium of a system with an equally large 



macroscopic subsystem of identical nature and dropping the subscripts and 
introducing the normalisation factor we get 



p = exp[-a-/3H-(y I N)H 2 ...] (27) 

P = 1 (28) 

kT 
-k/5 2 

y = — — (29) 
2G 

where C7 is the input specific heat at constant volume per particle The ensemble 
density is first expressed as a function of H and while integrating or summing over 
the appropriate states it is multiplied by g(E). Hence Eq. (27) is expressed as a 
function of H though the earlier equations in this section are all a function of E, 
referring as they do to microcanonical calculation. Eq. (27) represents the 
microcanonical ensemble and it can be shown from higher derivatives of Sr that 
every term in the exponent is linear in N. No other definition of equilibrium i.e. 
Ns/Nr^I would have lead to the microcanonical density for the system. For a simple 
system like that of N harmonic oscillators where the density of states is known one 
could verify that Eq.(27) indeed represents the reciprocal of density of states. By 
specifying the first few LaGrange multipliers we are in effect trying to approximate 
l/g(E). The canonical ensemble has one input parameter T, while the inclusion of 
second order term in H calls for another input parameter, C/ . The whole purpose of 
this exercise is to show through the intermediate ensemble 

p = exp[-a-f3H-(y/ N) H 2 ] po) 

that the average thermodynamic properties obtained from all these ensembles need 
not be the same even in the thermodynamic limit except at high temperatures. The 
justification or the lack of it for neglecting higher order terms in H would depend 
upon the temperature. From the form of the density it is easy to see that the higher 
order terms would have significant effect at low temperatures. In the next section we 
give a recipe for calculating the thermodynamic properties from all these ensembles. 



DEVALUATION OF THERMODYNAMIC QUANTITIES 

Define the new partition function for the intermediate ensemble as 

Z = Qxp(a)=\dr exp[-/3H-(y/N)H 2 ] qd 

Then by Eq. (7) 

S P = k[a + p(H) + (/ / Af)(# 2 )] (32) 

From the above equations, replacing the average energy of the ensemble by U we get 

- kT In Z = U - T\S P - k{y I N)(h 2 (33) 

For any other intermediate ensemble we have 



m 



-JcT]nZ = U-T[S P -kL filH 1 j] (34) 

Comparing Eq. (33) and Eq. (34) with the thermodynamic relation 

F = U-TS (35) 

we define the free energy and the thermodynamic entropy as 

F == -kT In Z (36) 

S = S P -l^(H^ (37) 

the average energy of the ensemble being the natural choice for average 
thermodynamic energy. That exp(-fiF) always represents the partition function could 
be easily interpreted as the maximisation of volume in phase space representing the 
minimisation of free energy in an appropriate way. Eq.(37) represents the only 
choice of 5 that would lead to thermodynamic equilibrium between two systems 
when 57+52 is maximised. This will be discussed in section VI. For any p choose 5 
such that this condition holds good. Thus in Eq.(36) and Eq.(37) no other regrouping 
of terms could be chosen to represent the thermodynamic entropy. 

The other thermodynamic quantities cannot in general be expressed as 
derivatives of free energy. The reason lies in the following. The thermodynamic 
relation 

n 

TdS = dE+Y, Xi dxi (38) 
(=i 

where Xi's including E and 5 are the state variables and the partial derivatives Xi's 
are the conjugate forces is a micro law valid for every member in the ensemble with a 
given energy E. It arises from 

T[(S+dS) -S] = [(E+dE) -E]+i X{(x,+dx,)-x,] (39) 

(=i 

Using Ehrenfest principles 6 if we want to calculate Xi from the relation 

Xi = \_SE I SXi\(x\, X2,.. xa - i), x (i + i).... x„), S (40) 

where the constancy of 5 refers to the constancy of the density of states 
at a given energy E, we have to bear in mind that systems with a given energy E 
occur with a definite probability. Hence, on multiplying Eq. (40) by p and integrating 
between appropriate limits, we get, in the energy representation 

oo 

X) = exp(— a)\dT— [dE / dxi](», X 2, ..x„ + g(E)p(E) (4i) 



By multiplying Eq. (38) by the probability density and integrating between 
appropriate limits we obtain 

dE) = (TdS)-L{Xdx l ) (42) 



from which it is not possible to derive Eq. (41) and hence there is no way one could 
derive 



d(E} = (d(Eyd(S})v.»rf(sy(d(EydV^ 

(43) 

which is what is used in deriving the conjugate forces in the canonical ensemble 
formulation. It means that although in Eq. (39) the infinitesimals can be taken about 
the average values of state variables, with 5 and E replaced by their average values, 
the conjugate forces, which are given by the partial derivatives of the former have to 
be derived only from Eq. (41). So (Xi) need not in general satisfy the relation 

Xi = [d U I dxi] («,... s) (44) 

though it may be valid for certain cases like the canonical ensemble. For example in 
the case of the intermediate ensemble pressure is given by 

P = r o dE-[dE/dVl NS) g(E)exp(-a-pE-(y/N)E 2 ) 

(45) 

and is not equal to 

d/dV[-kT\nZ] (46) 

Eq. (45) alone ensures the validity of virial theorem: PV=2/3U for free particles 
which can be derived from mechanical considerations alone and Eq. (44) violates it. 
In short the processes of ensemble averaging and Legendre transformation need not 
commute. Ensemble average has to follow Legendre transformation and not precede 
it. The recipe hence is to apply Eq. (38) to every microcanonical ensemble and use 
Eq. (36), Eq. (37) and Eq. (41) for evaluating the thermodynamic properties. 

To complete the analogy between thermodynamics and statistical mechanics 
one has to show that the thermodynamic limit exists. With the intermediate ensemble 
for a hard core finite range attractive potential the existence of the thermodynamic 
limit can be proved along the lines discussed in Thopmson 7 for Gibbs ensemble. In 
all these derivations the laws of thermodynamics are kept in tact. The transition from 
statistical mechanics to thermodynamics is made through the axiom given by Eq. (21) 
and whatever we derive with different densities has to be and is consistent with 
thermodynamics. These calculations therefore do not represent a rederivation of the 
macroscopic thermodynamic laws from microscopic mechanics. 

V.AVERAGE ENERGY AND FLUCTUATIONS 

In the energy representation, partition function for Eq.(30) can be represented by 

oo 

Z = jdE g(E) exp[-/3E-(y/N)E 2 ] m 





Taylor expansion of lng(E) yields, 



Z = g(E)UdEexp{-l3(E-E)-(y/N)(E-Ey- 



i=i i\ dE' 

(48) 

It is obvious that the integrand does not always have a saddle point as it depends on 
the signature of the coefficients. In evaluating the canonical partition function 

oo 

Z = \dE g{E) exp(-j3£) (49) 





(50) 



is claimed to be the saddle point equation. It is not the saddle point equation but 
the defining equation for E(T) in the microcanonical ensemble calculation. 
Hence 

E-E-Tk In g(E) (5i) 

is the relation between averages obtained from microcanonical ensemble. It is 
important to note that F is not given by the partition function in by Eq. (49). 
Comparing the partition functions given by Eq. (48) and Eq. (49) we infer readily that 
for a given T average energy in the intermediate ensemble is smaller than that in the 
canonical ensemble and as we conserve higher and higher powers of H in the 
ensemble density, the value of E decreases steadily and we get the lowest value in 
the microcanonical case. From the nature of the higher order multipliers which are all 
T dependent, we also realise that the change in E(T) is not one of mere scaling or 
shifting. Since the LaGrange multipliers are all inversely proportional to T, we 
conclude that when at high temperatures all these ensembles must yield the same 
results they may not yield the same results at low temperatures even in the 
thermodynamic limit, for, all the terms in the exponent area linear in N. This analysis 
does not exclude the possibility of certain systems giving the same results in 
microcanonical and canonical ensembles, like for instance classical ideal gas. In this 
case the density of states obeys a power law and we also know that both the canonical 
and microcanonical ensembles yield the same results. In such a case as this, the 
limiting ensembles yield the same results because the saddle point exists in the 
partition function. That is, by solving the canonical partition function by the saddle 
method we are in effect doing a microcanonical ensemble calculation. In other cases 
through the Hamiltonian and LaGrange multipliers we are in effect trying to 
approximate the density of states. In some cases there may not be a saddle point in 
the partition function but the average energy might have as in 



where c is a constant, though this might not represent any realistic case. Even in 
those cases where there is no saddle in the expression for average energy, the 
statistical entropy can always be solved by saddle point approximation for a suitable 
upper limit m in the truncated Eq. (13). In this case a microcanonical calculation 
would lead to a negative temperature as Eq. (50) does. But the canonical ensemble 
yields a relation anlogous to Eq. (35) ! 

That the average energy, entropy and free energy are related by a Legendre 
transformation - except for a difference due to terms of order InN negligible in the 
thermodynamic limit is usually interpreted as the equivalence 8 of ensembles in the 
thermodynamic limit. We have shown here that given any density of states we could 
always find these thermodynamic quantities satisfying such a relation, from the axiom 
given by Eq. (21) and that these averages Xi(T) vary from ensemble to ensemble. 
Thus this relation is necessary but not sufficient. As the densities differ, the partition 
function, and hence the average energy and free energy differ and hence does 
thermodynamic entropy and the difference is more than in terms of order InN. 

E~ 2 -E 2 

It is easy to show that goes inversely as ./V in all these ensembles. The 

E 2 

second central moment decreases steadily as we move from canonical to 
microcanonical ensemble. In the case of the intermediate ensemble we have 

((£-£) 2 ) = {(E-E c ) 2 y(E-E c y (53) 

where E c is the average energy in the canonical ensemble and within the brackets 
representing the averaging with respect to the intermediate ensemble it remains a 
constant. Hence the second central moment in the intermediate ensemble is smaller 
than that in the canonical ensemble. The width decreases with the inclusion of higher 
powers of H. The specific heat per particle C obtained from the intermediate 
ensemble includes correction terms besides C/. It is given by 



kT 2 NC R = (H 2 y(H) 2 - kT 2 d I dT[y I N {(h 3 } - (H)(h 2 } } ] 

(54) 

The grand canonical ensemble yields the same equation for (Xi(T)/N) as the original 
ensemble. But if we want to obtain (N(T)) it is important to include, for reasons 
similar to those mentioned with regard to the variable E, higher order terms in ./V and 
also cross terms involving E and N. 



VI. COUPLING BETWEEN SYSTEMS 

Consider two systems of energy Ej and E 2 and particle number Ni and N2 
respectively. If the two systems are coupled to form a composite system isolated 
from the surroundings, then the composite system would be in equilibrium and the 
sum of the statistical entropies is a maximum. In general for any variation in the 
energies of a composite system of n sub-systems 



s[j:s p ]=o (55) 

i=\ 

Since the total energy is fixed 

£<5 Ei=0 (56) 
(=i 

Condition in Eq. (55) corresponds to allowing maximum number of states for the 
composite system. If we use the intermediate ensemble to describe the systems, using 
Eq. (32) and Eq. (56) we could reduce this condition for the case n=2 to 

Q3, - pi) + 2[(/i / Ni)Ei - (Y2 1 N2)E 2 ] = (57) 

This condition for statistical equilibrium does not imply thermodynamic equilibrium. 
If the two systems are of the same nature, then, at thermodynamic equilibrium, it 
further demands that 

Ei Ei 

= (58) 

Ni N 2 

This, nothing but the extensive property of energy, will be true by very definition of 
equilibrium. If the two systems are described by microcanonical or canonical 
ensemble the condition for statistical equilibrium is the same as that for 
thermodynamic equilibrium. But if the composite system is made up of three or more 
systems none of these ensembles stipulates thermodynamic equilibrium as the 
condition for statistical equilibrium as can be readily inferred from Eq. (55) and Eq. 
(56). These would be valid in general because of the maximisation of S p and the 
nature of the condition would depend on what ensemble we use to describe the 
system. It is easy to verify that statistical equilibrium is a transitive property. 
Physically Eq. (54) means that if two systems at different temperatures were brought 
together the composite system when shielded from the surroundings would be in 
statistical equilibrium even before thermodynamic equilibrium is attained. Thus 
thermodynamic equilibrium is a special case of statistical equilibrium. 

Finally at statistical equilibrium we should be able to cast the density in the 
same form as p. This means that we have to find an effective, and not equilibrium f$ 
and y/N in terms of and ji. Starting with the expression for total statistical 

entropy and using the relation ( E? ) — (Ei) on the assumption that the integral 



for the saddle point could be solved by saddle point approximation, we get 



S P = S PI + S P2 = k[cCi+frEi + (yJ Ni)Ei 2 +a2+P2E 2 +(Y2/ N 2 )E 2 2 ] 

(59) 

and using the condition for statistical equilibrium given in Eq. (55) we cast S p into the 
same form as S P , 

S P =k[a+kfr+foE+\{^+^W] (60) 
2_ 4 N, N, 

where we have neglected (E i — E 2) 2 in comparison with (£" 1 + E 2) 2 . The 
effective values of f$ and ( /A 7 ) are given by the coefficients in Eq. (60) and 

p = cxp[-a- t ^H-{^+-^- }H>] (6i) 

is a stationary state of Liouville's equation. Eq. (61) represents the state of the 
system at t=to when the individual systems are brought together and is the initial state 
of the composite system whose time evolution is dictated by Liouville's equation. 
The final equilibrium state of the composite system will also be given by Eq. (60) 
except that fa and fa would both correspond to the temperatures at which the two 
systems equilibrate. 



VII. CONCLUSION 

Through a comparative study of ensembles ranging from canonical to 
microcanonical, we have proved that all these ensembles may not yield the same 
results even in the thermodynamic limit except at high temperatures. While a simple 
system like that of a system of harmonic oscillators may be investigated to verify 
these results and check the qualitative and quantitative differences, it must be pointed 
out that it is only an academic exercise. The more interesting cases are those where 
we do not know the density of states and try to approximate it through experimentally 
measured parameters. As a byproduct of this analysis we have obtained a general 
form for probability density function in an interval. Power series expansion for 
functions of several variables is derived and convergence questions addressed using 
spherical polar coordinates. 
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APPENDIX 

A probability density of the form given in Eq. (19) can be used when solutions to 
differential equations are sought for. Before deriving it in yet another way power 
series in several real variables is considered below. It is well known 9 that starting 
from 

Pn(x) = Y,lo ajC ' (62) 

where 

a = 1 / i ! [d 1 1 dx l P(x)] ( , = 0) (63) 

Taylor and his pupils defined power series starting from 



fix) = LIo*' li] \- di ldx l f{x)l* = 0) + R N (64) 
by taking the upper limit to infinity. Analogously if one begins from 10 

and defines 

aj = -^—[d ii+j) Idxidy'P]^ = 0> , = 0) (66) 

one could show for a and A: real, since 

(l+ax) k = 1+ k(ax)/ 11+ k(k-l)(ax) 2 / 21+ (67) 

is convergent for laxl < 1, 

(1 + xy y = l + k(xy) l\\+k(k- l)(xy) 2 / 2 !+ (68) 

is convergent for \xy\<l or equivalently r<l in polar coordinates, where (x,y)eR 2 . 
Similarly since the power series for exp(a x) is convergent for all xeR. 

exp(xy) = 1+ (xy) / 1!+ (xy) 2 / 2!+ +{xy) r I r\+.... (69) 

and the above series is convergent for all (x,y)eR 2 . Taylor expansion for two 
variables is defined by 

fi^^l^a^y^R, (70) 

and the coefficients in (68) and (69) and (70) are precisely those derived from 

an = —[d Idxtdytf]* = o,, = 0) (7i) 

Hence power series in two variables can be defined as 

f(x,y)=LZ ) , m a>x'y' 

if convergence questions for functions of several variables as in Eq.(72) are 
addressed using polar coordinates. Equation (72) is the limiting form of Eq.(70) 

and Eq. (71) gives a more natural expression for a t j 's. Similarly the power series 
expression for logarithm could also be derived to arrive at Eq. (19). Extension to 
several variables is straightforward. Extrema of functions of several variables could 
also be found easily solving just one equation in the radial coordinate and checking 
the signature of the second derivative in that coordinate 11 . 



